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Abstract 

We study the problem of extracting information about composition and equation of state of dense 
matter in neutron star interior using axial w-modes. We determine complex frequencies of axial 
w-modes for a set of equations of state involving hyperons as well as Bose-Einstein condensates of 
antikaons adopting the continued fraction method. Hyperons and antikaon condensates result in 
softer equations of state leading to higher frequencies of first axial w-modes than that of nuclear 
matter case, whereas the opposite happens in case of damping times. The presence of condensates 
may lead to the appearance of a new stable branch of superdense stars beyond the neutron star 
branch called the third family. The existence of same mass compact stars in both branches are 
known as neutron star twins. Further investigation of twins reveal that first axial w-mode frequen- 
cies of superdense stars in the third family are higher than those of the corresponding twins in the 
neutron star branch. 

PACS numbers: 97.60.Jd, 26.60.-c, 04.40.Dg 
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I. INTRODUCTION 



Shortly after the discovery of pulsars, the study of dense matter in the core of neutron 
stars gained momentum With the advent of satellite based observatories such as Einstein, 
ROSAT, Hubble space telescope XMM-Newton, and Chandra x-ray observatory, the study 
of neutron stars has entered into a new era. Observations using these facilities as well as 
other facilities are pouring in very exciting data on neutron stars. It is now possible to 
estimate masses, radii, moment of inertia and surface temperatures of neutron stars from 
the observations. Those findings ,in turn, may shed light on the gross properties of cold 
dense matter far off from normal nuclear matter density in neutron star interior. Such a 
cold dense matter in the core of a neutron star, so far, has not been produced in laboratories 
on the earth. Therefore neutron stars are very useful laboratories for the study of highly 
dense matter in its core. 

On the other hand, there is a growing interplay between the physics of dense matter in 
relativistic heavy ion collisions and that of neutron stars. Though the Quantum Chromo- 
dynamics predicts a very rich phase structure of dense matter, we can only probe a small 
region of it in the laboratories. The study of dense matter in heavy ion collisions reveals 
many new and interesting results such as the modifications of hadron properties in dense 
medium, the properties of strange matter including hyperons and (anti)kaons, the strength 
of attractive antikaon-nucleon interaction and the formation of quark-gluon plasma. These 
empirical informations from heavy ion collisions may be useful in understanding dense mat- 
ter in neutron star interior. Heavy ion experiments in the upcoming Facility for Antiproton 
and Ion Research (FAIR) at GSI, Germany would produce dense matter with baryon density 
a few times normal nuclear matter density and temperature a few tens of MeV and provide 
us important information about dense matter. 

Gravitational wave astrophysics opens a new window to probe the neutron star interior 
0]. Neutron stars are accompanied by strong gravitational field. Non-radial oscillations of 
neutron stars generate gravitational waves. Such a situation may arise in neutron star glitch, 
merger of compact stars, a phase transition in the core and gravitational collapse of a massive 
star. Some large interferometric gravitational wave detectors such as LIGO in USA, VIRGO 
in Italy and GEO600 in Germany are now operational. A neutron star has a large number 
of families of pulsation modes [3, 3). Those modes are classified according to restoring 
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forces acting on a perturbed fluid element. Important modes among them are fundamental 
f-mode associated with global oscillations of the fluid, g-mode due to buoyancy and p-mode 
due to pressure gradient. General relativity also predicts w- modes which are space-time 
modes. Unlike other pulsation modes of neutron stars, w-modes do not have counterparts in 
Newtonian gravity. Equations governing perturbation of a static and non-rotating neutron 
star are decomposed into two sets depending on how perturbations transform under parity. 
Polar perturbation corresponds to even parity ( — 1)' and axial perturbation relates to odd 
parity (— . Various groups (5], [f], Q, ^, ijl investigated w-modes in great details and 
found interesting varieties of axial w-modes |7|, |8|]. It is to be noted that w-modes have 
higher frequencies than those of /, g and p modes. 

In general relativity, oscillations of neutron stars are damped by gravitational radia- 
tion. Consequently, the corresponding oscillatory modes are quasi-normal modes and their 
frequencies are complex. Once gravitational waves from a neutron star are detected, the 
spectral analysis of the detected signal would tell us the frequency and damping time of 
the quasi-normal mode responsible for gravitational waves. Now the question is how we 



determine the mass, radius and properties of dense matter from the frequency anc 
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damping 



time of the mode. This is achieved by solving the "inverse problem" [3], [k 
solve this problem, we have to compute frequencies and damping times of modes that could 
generate strong gravitational waves using different compositions and equations of state and 



compare them with those obtained from observations. This way we could probe the neutron 
star interior. 

Various novel phases of matter might appear in neutron star interior because of large 
baryon and lepton chemical potentials there. Glendenning jl] predicted the presence of 
hyperons in neutron stars. Kaplan and Nelson discussed Bose-Einstein condensation of K~ 
mesons in neutron star matter due to strongly attractive antikaon-nucleon interaction 13]. 
The formation of quark matter is another possibility in neutron star core In this paper, 
we are motivated to investigate the role of hyperons and antikaon condensation on frequencies 
of quasi-normal modes of neutron stars. In particular, we would investigate how the strength 
of antikaon optical potential in nuclear matter which is crucial for antikaon condensation in 
neutron stars, affects axial w-modes. In this connection, we want to calculate the frequency 
and damping time of axial w-mode of a neutron star including Bose-Einstein condensate of 
negatively charged kaons for a set of values of antikaon optical potential depth. Further we 
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investigate the role of hyperons on frequency and damping time of axial w-modes. 

This paper is organised in the following way. The formalism is discussed in section II. 
Results are explained in section III. Section IV provides summary and conclusions. 



II. FORMALISM 

Chandrasekhar and Ferrari showed how the problem non-radial oscillations of a compact 
star could be reduced to the scattering of incident gravitational wave by the static space- 

' neutron 
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time of spherically symmetric star [5|, [6[ . Here we are interested in axial w-modes o: 
stars. The axial perturbation of the compact star is described by the metric 



ds 2 = e 2v (dtf - e 2i '(d<j) - q 2 dr - q 3 d9 - udtf - e 2 ^(drf - e 2 ^(d6) 2 , (1) 

where metric functions u, if), ^ and ^ retain their unperturbed values and to, q 2 and 
give rise to the perturbations. The radial evolution of equations for axial perturbations is 
described by the following by the Schrodinger like equation with a potential barrier V(r) 

^ + u) 2 -V(r))z l = 0, (2) 



dr 



e 2v 



V(r) = — (l(l + l)r + r 3 [e(r)-P(r)]-QM(r)) , (3) 

dv 1 dP 

dr ~ ~(e(r) + P(r))~d~r~ ' ^ ' 

where the tortoise coordinate is r* = J r e-'+^dr, M(r) is the mass enclosed in a radius 
r, uj is the complex frequency of the mode and e and P are energy density and pressure 
respectively. Both energy and pressure outside the star become zero and the Eq. (j2J) reduces 
to the Regge- Wheeler equation with 

™=(l-^)(^-£). ( 8) 

and r* = r + 1Mln(r — 2M). We determine frequencies of quasi-normal modes by solving 
Eq. (j2J) with appropriate V(r) inside and outside the neutron star and boundary conditions 
such that the solution is regular at the center, continuous at the surface and behaves as a 
purely outgoing wave at infinity. 

Now we focus on the construction of equation of state (EoS) (pressure versus energy 
density) within relativistic field theoretical models. In this paper, we consider EoS with 



different compositions and values of incompressibility of nuclear matter. Firstly we discuss 
the EoS including neutrons (n), protons (p), hyperons (H=A,£ + ,X _ ,X ,H~,H ) and electrons 
(e) and muons(/z). The baryon-baryon interaction is mediated by the exchange of a, u), p 
mesons and for hyperon-hyperon interaction, two additional mesons - scalar meson / (975) 
(denoted hereafter as a*) and vector meson 0(1020) [15) are incorporated. Therefore the 
Lagrangian density for baryons in the charge neutral and /9-equilibrated hadronic phase is 
given by 

Cb = ^^B («7m^ ~ m B + 9aBCT - g^Bl^ - QpBl^B ■ p'') ipB 
B 

--u^ + \mluj^ - ~ PfX1/ ■ pT + Ur? pPil ■ p» + C YY . (6) 

The isospin multiplets for baryons B = N, A, £ and H are represented by the Dirac spinor 
ty B with vacuum baryon mass mg and isospin operator t B and u^ u and p^ v are field strength 



tensors. The scalar self-interaction term 



ld| is 



U(a) = l -g 2 o* + \g^ . (7) 
The Lagrangian density for hyperon-hyperon interaction (Cyy) is given by 
Cyy = J2^ b (9**bV* - 94>Bl^) ^b 

B 

+~ (d»a*d»a* - mW) - \^ V <T + \™\^ ■ (8) 



As nucleons do not couple with strange mesons, g a *N = g<j>N = 0. 

Next we discuss equations of state undergoing first order phase transitions from nuclear 
to K~ condensed matter denoted by npX~ and hyperon matter to K (K~, K°) condensed 
matter denoted by npHAf. Here we have pure hadronic phase described by Eq. (jSJ) and 
antikaon condensed phases and the mixed phase of two pure phases. The constituents of 
/3-equilibrated matter in both phases are neutrons, protons, hyperons, electrons and muons. 
Baryons are embedded in the condensate in the condensed phase. We describe the interaction 
of (anti)kaons in a relativistic field theoretical model. The Lagrangian density for (anti)kaons 



in the minimal coupling scheme is 
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2l| 



C K = D*KD»K - mfKK , (9) 



where the covariant derivative is = + ig^K^^ + ig^K^n + WpK^K • P M and the effective 
mass of (anti)kaons is m* K = rriK — Qokv — 9o*k<3* ■ The mixed phase of hadronic and K~ 
condensed matter is governed by the Gibbs conditions for thermodynamic equilibrium along 



with global charge and baryon number conservation laws [22]. Working in the mean field 
approximation, we obtain energy density and pressure in each pure phase and those are 
given by Ref. 21]. 



III. RESULTS AND DISCUSSION 



Nucleon-meson coupling constants obtained by reproducing saturation properties of nu- 
clear matter for two values of incompressibility K = 240 and K = 300 MeV, are taken 



from Ref. 



231 ] . These are known as GM sets. Hyperon- vector meson coupling constants 
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24 



251 ] . The scalar meson a 



are determined from SU(6) symmetry of the quark model 
meson coupling to hyperons is calculated from the potential depth of a hyperon (Y) 



UP (n ) 



(10) 



in normal nuclear matter. Potential depths of hyperons in normal nuclear matter are ob 
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tained from the analysis of hypernuclei data 2 



271 . |28j| . Similarly hyperon-a* coupling 
constants are estimated by fitting them to a potential depth, \no), for a hyperon (Y) 
in hyperon matter. We adopt the values of hyperon potential depths for this calculation 
as quoted in Ref. 21] . Next we determine kaon- vector meson coupling constants using the 



quark model and isospin counting rule JjJ • The scalar coupling with kaons is obtained from 



the real part of K optical potential depth at normal nuclear matter density [29|, |3 



(11) 



Ur K) 



where a and loq are mean values of the meson fields. We perform our calculation for an- 
tikaon condensation with K = 300 MeV and two values of antikaon optical potential depth 
Ur^o) — —120,-160 MeV. Values of q n K are recorded in Ref. [21] . Coupling constants 



g a *K and g^x are obtained from Refs.[l5|, \21\\ . 



Now we present our results for different compositions and equations of state. Gravi- 
tational mass for static neutron star as well as superdense star sequences are shown with 
radius in Figure 1 for different values of incompressibility. For np matter we find the EoS 
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with K = 300 MeV is stiffer than that of the case with K = 240 MeV. Consequently the 
maximum mass is larger in the former case. With the appearance of hyperons and/or an- 
tikaon condensates, equations of state become softer. For npA matter and K = 300 MeV, 
the threshold density of A hyperons is 2.3t7, . On the other hand, for npK~ matter with 
K = 300 MeV, K~ condensation sets in at 2.23no and the mixed phase terminates at 3.59n 
for Ux(no) = -160 whereas the phase transition for U^(no) = -120 MeV begins at 3.05rio 
and there is no mixed phase in this case. Stronger the attactive antikaon optical potential 
depth, softer is the EoS. Softer EoS involving hyperons and antikaon condensates lead to 
reduction in maximum masses than that of the np matter as evident from Fig. 1. For 
npHK~ K° case with K = 300 and Uj^(n ) = —160 MeV, A hyperons appear at 2.51r2, 
after the onset of K~ condensation at 2.23n . The mixed phase ends at 4.0n in this case. 
Further K° condensation which is treated here as a second order phase transition, begins at 



4.06™ 



211 ]. Heavier strange baryons are populated beyond 6no 2lj. For npHK K° case, it 



was noted that after the stable neutron star branch, there was an unstable region followed 
by another stable branch of compact stars 2lJ. As white dwarfs and neutron stars form 
first two families of compact stars, this branch of superdense stars beyond the neutron star 
branch is called the third family of compact stars [3jJ. Further the existence of non-identical 
stars having same mass could be possible because of partial overlapping of the neutron star 
and third family branches. These pairs of compact stars are known as "neutron star twins" 



32 



33j . Superdense stars in the third family branch have different compositions and smaller 



radii than their counterparts in the neutron star branch 21]. 

Next we calculate the frequency (v) and damping time (r) from Eq. (j2J) using equations 
of state described above. We exploit both the continued fraction method jlOj and integration 
method [7] to find complex eigenfrequency of the Schrodinger like equation (T5]). Here we 
present our results obtained by using the continued fraction method. The frequency and 
damping time are related to the real and imaginary parts of eigenfrequency through v = 
^(Muj ) kHz and r = 4.937 x //s, where n=M/M and u , w u M and M are 
measured in units of km. Figures 2 and 3 display frequency and damping time of first axial 
w-mode as a function of neutron star compactness (M/R) respectively. In both figures, the 
dashed line corresponds to np matter with K = 240 MeV whereas the solid line represents 
np case with K = 300 MeV. We note that the softer EoS leads to higher frequency of first 
axial w-mode for each M/R value as shown in Fig. 2 whereas it is the reverse for damping 
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time as evident from Fig. 3. 

In the next paragraphs we discuss the role of exotic matter on frequency and damping 
time of axial w-mode. Firstly we study axial w-modes of neutron stars involving npA matter 
along with those of np matter for K = 300 MeV. Frequencies and damping times of first 
axial w-mode for both cases are exhibited in Figures 4 and 5. The appearance of A hyperons 
makes the EoS softer compared with the np case. This leads to higher frequencies of first 
axial modes for oscillating neutron stars including A hyperons as shown by the upper curve 
in Fig. 4 whereas damping times in this case are given by the lower curve in Fig. 5. 
Frequencies and damping times of first three I = 2 axial w-modes corresponding to the EoS 
involving npA matter with increasing compactness are recorded in Table I. 

We continue our investigation of first axial w-modes of neutron stars including a first 
order K~ condensate. These results are shown in Figures 6 and 7. We have considered two 
values of antikaon optical potential depth at normal nuclear matter density = -120, -160 
MeV. For = -160 MeV, frequencies of first axial w-modes for neutron stars having a K~ 
condensate are shown by the top curve whereas those of npK~ EoS for Ux(n ) = -120 MeV 
have lower frequencies. This may be attributed to the fact that the EoS with (n ) = 
-160 is softer than that of Ux(n ) = -120 MeV. Similarly we find that damping times for 
Ur^o) = -160 MeV are reduced than those of U^(n ) = -120 MeV and np case. Frequencies 
and damping times of first three / = 2 axial w-modes of oscillating neutron stars including 
a K~ condensate with increasing compactness are recorded in Table II. 

Frequencies and damping times of first axial w-modes of neutron star as well as third 
family sequences as a function of compactness are shown in Figures 8 and 9. In Fig. 8 fre- 
quencies of first axial w-modes corresponding to the neutron star and third family branches 
are denoted by the bottom and top curves respectively. Further twins of superdense stars 
in the third family branch (top curve) are represented by filled squares in the neutron star 
branch (bottom curve). Keeping baryon number fixed, the collapse of an ordinary neutron 
star to its twin in the third family branch might release huge amount of energy. A part of it 
may be carried away by gravitational waves. This could be a good signal to probe interior 
of superdense stars. Damping times of superdense stars corresponding to the third family 
branch have higher values than those of the neutron star branch. 

We compare real and imaginary parts of eigenfrequencies of first axial w-modes for EoS 
including npA and EoS involving npK~~ with Ux(n ) = -160 MeV in Figures 10 and 11. 
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We find that values of real and imaginary frequencies corresponding to both EoS are close 
beyond compactness 0.18. We conclude from these results and those recorded in Table I and 
II that these equations of state are indistinguishable in axial w-modes. 

IV. SUMMARY AND CONCLUSIONS 

We have studied quasi-normal modes of oscillating neutron stars with different com- 
positions and equations of state. This problem has been investigated adopting continued 
fraction method. We have computed frequencies and damping time scales of axial w-modes 
corresponding to EoS involving np, npA and npK~ and npRK'K matter. We find that 
hyperons and/or antikaon condensates result softer EoS which ,in turn, lead to higher fre- 
quencies and lower damping time compared with those of np case. Axial w-modes might 
distinguish between exotic matter and nucleons-only matter. However, it would be difficult 
to differentiate the EoS of hyperon matter from that of antikaon condensed matter using 
axial w-modes as probe. Further our study has revealed that freuqencies of superdense stars 
in third family branch are higher than their twins in the neutron star branch. 
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TABLE I: First three values of characterisitic 1=2 axial w-mode frequencies and damping timescales 
for neutron stars made of n, p, A and electrons and muons. 



M R § v t 
(M ) (km) (kHz) Qus) 



1.822 12.892 0.2 6.830 79.571 
10.055 47.711 
13.391 50.397 



1.867 12.325 0.22 6.796 87.129 
9.8556 52.440 
12.800 51.996 



1.885 11.514 0.24 6.800 94.766 
9.595 60.417 
12.148 54.029 
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TABLE II: The first three values of characteristic 1=2 axial w-mode frequencies and damping 
timescales for neutron stars undergoing a first order phase transition from nuclear to K~ condensed 
matter for antikaon optical potential depth at normal nuclear matter density U^(no) = — WOMeV. 



M R § v t 
(M ) (km) (kHz) (fis) 



1.680 12.401 0.2 6.996 69.554 
10.107 46.401 
13.691 51.040 



1.770 11.641 0.22 6.941 80.302 
9.832 52.712 
12.779 53.119 



1.825 11.068 0.24 6.900 89.476 
9.593 60.270 
12.149 54.333 
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Fig. 1. Gravitational mass is plotted with radius for compact stars having different compositions 
as explained in text. 
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Fig. 2. Frequency of first axial w-mode is plotted as a function of neutron star compactness for 
nucleons only matter corresponding to incompressibility K=240 and 300 MeV. 
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Fig. 3. Damping time of first axial w-mode is plotted as a function of neutron star compactness 
for nucleons only matter corresponding to incompressibility K=240 and 300 MeV. 
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Fig. 4. Frequency of first axial w-mode is plotted as a function of neutron star compactness for 
np as well as npA matter with K=300 MeV. 
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Fig. 5. Damping time of first axial w-mode is plotted as a function of neutron star compactness 
for np as well as npA matter with K=300 MeV. 
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Fig. 6. Frequency of first axial w-mode is plotted as a function of neutron star compactness for 
np and npK~ matter with K=300 MeV and U R {n ) = -120, -160 MeV. 
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Fig. 7. Damping time of first axial w-mode is plotted as a function of neutron star compactness 
for np and npR- matter with K=300 MeV and U R (n ) = -120, -160 MeV. 
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Fig. 8. Frequency of first axial w-mode is plotted as a function of compactness for neutron stars 
including n,p,A,E-,E°, Yr ,K~ and K° with K=300 MeV and U R (n ) = -160 MeV. The lower 
curve shows the neutron star branch whereas the third family branch is shown by the upper curve. 
Neutron star twins corresponding to the superdense stars in the third family branch are denoted 
by filled squares in the lower curve. 
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Fig. 9. Damping time of first axial w-mode is plotted as a function of compactness for neutron 
stars including n,p,A,E~,E°, Yr ,K~ and K° with K=300 MeV and U R (n ) = -160 MeV. The 
lower curve shows the neutron star branch whereas the third family branch is shown by the upper 
curve. Neutron star twins corresponding to the superdense stars in the third family branch are 
denoted by filled squares in the lower curve. 
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Fig. 10. Real part of eigenfrequency of first axial w-mode is plotted as a function of compactness 
for npA and npR- matter with K=300 MeV and U R (n ) = -160 MeV. 
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Fig. 11. Imaginary part of eigenfrequency of first axial w-mode is plotted as a function of 
compactness for npA and npK~ matter with K=300 MeV and U^(no) = -160 MeV. 
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